Testing Bell's inequality in constantly coupled Josephson circuits by effective single- qubit operations 
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In superconducting circuits with interbit untunable (e.g., capacitive) couplings, ideal local quantum opera- 
tions cannot be exactly performed on individual Josephson qubits. Here we propose an effective dynamical 
decoupling approach to overcome the "fixed-interaction" difficulty for effectively implementing elemental logi- 
cal gates for quantum computation. The proposed single-qubit operations and local measurements should allow 
testing Bell's inequality with a pair of capacitively-coupled Josephson qubits. This provides a powerful ap- 
proach, besides spectral-analysis [Nature 421, 823 (2003); Science 300, 1548 (2003)], to verify the existence of 
macroscopic quantum entanglement between two fixed-coupling Josephson qubits. 



PACS number(s): 03.65.Ud, 03.67.Lx, 85.25.Dq. 



I. INTRODUCTION 



Non-locality (i.e., entanglement) is one of the most pro- 
found features of quantum theory and plays an important role 
in quantum information processing, including quantum com- 
putation and quantum communication |JJ. Mathematically, an 
entangled pure state of a composite system is defined as a state 
that cannot be factorized into a direct tensor product of the 
states associated with individual subsystems. The presence 
of multi-partite entanglement is necessary for implementing 
quantum algorithms that are exponentially faster than classical 
ones 1 1]. Also, entanglement may offer some new informa- 
tion transfer modes, e.g., teleportation and quantum cryptog- 
raphy 1 1 ] . Therefore, generating and verifying entanglement 
between qubits are of great practical importance |2]. 

The existence of entanglement can be verified by using var- 
ious experimental methods, e.g., quantum tomographic tech- 
niques. Bell-state analysis, quantumjump measurements, etc. 
(see, e.g., |3, 4]). Also, spectral analysis has been used to 
probe the existence of two-qubit entanglement in coupled 
Josephson qubits However, the degree of entanglement 

between the two always-present interacting qubits changes 
rapidly |7] and, at certain times, two-qubit states can be al- 
most separable. Thus, verifying the instantaneously generated 
entangled state in coupled systems |8], and using it to realize 
quantum information processing, e.g., teleportation and quan- 
tum memory, are very important challenges. 

Historically, Bell's inequality always served as one of the 
most important witnesses of entanglement: its violation im- 
plies that entanglement must be shared by the separate parts. 
Numerous experimental tests of Bell's inequality have been 
made with entangled photons separated far apart (e.g., up to 
500 m) 1 9] and entangled closely-spaced trapped ions (e.g., 
separated a few micrometers apart) \ 10]. A Bell-like inequal- 
ity has also been tested via single-neutron interferometry by 
measuring the correlations between two entangled degrees of 
freedom (comprising spatial and spin components) of single 
neutrons |11]. The results from these experiments strongly 
violate the tested Bell's inequalities, and thus agree with quan- 



tum mechanical predictions. Very recently, preliminary pro- 
posals have been explored for testing Bell's inequalities with 
switchable Josephson qubits 1 12|]. 

Almost all proposals for manipulating quantum information 
rely on the execution of both single-qubit and two-qubit gates. 
In some cases (e.g., trapped ions, QED cavities, and quantum 
dots O), single-qubit operations are easy to realize by ap- 
plying certain controllable local fields. However, the interbit 
couplings are fixed andjtncontrollable in current experimental 
Josephson circuits |4, 5, 6] and NMR systems |14]. Never- 
theless, qubits in NMR are individually addressable, because 
the coupling constants Jij between the ith and jth qubits 
are sufficiently weak (i.e., much smaller than the differences 
Awij = \uji — LOj I between the eigenfrequencies of the qubits, 
e.g., Jij / lS.ijJij < 10~'^ O). However, the usual capacitive 



coupling in Josephson circuits \4,'y\ is relatively strong, thus 
making it difficult to perform local single-qubit operations on 
individual qubits. The method fl3l of physically arranging 
qubits (i.e., using two or three physical qubits to encode a 
logical qubit) cannot be directly used in the present two-qubit 
system. 

Here we develop a new approach to overcome the seri- 
ous "fixed-interaction" difficulty and effectively implement 
desired single-qubit operations on selected qubits. This dy- 
namical decoupling method can be used to generate long-lived 
entangled states in coupled Josephson qubits. The long-lived 
entanglement obtained here, assisted by the proposed local 
single-qubit operations, should allow to test Bell's inequal- 
ity with a pair of capacitively coupled Josephson qubits. Its 
violation would provide another robust physical evidence of 
statistically nonlocal correlations at the macroscopic scale. 



n. DYNAMICAL DECOUPLING 

We consider the nano-circuit sketched in Fig. 1. This 
is similar to that in Ref. |5], just replacing the left Joseph- 
son junction there by a superconducting quantum interfer- 
ence device (SQUID) loop with controllable Josephson en- 
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FIG. 1 : Two capacitively coupled charge qubits. The quantum states 
of two Cooper-pair boxes (i.e., qubits) are manipulated by control- 
ling the applied gate voltages Vg^ , Vg^ and external magnetic fluxes 
$1, <1>2 (penetrating the SQUID loops). Pi and P2 (dashed line 
parts) read out the final qubit states. 



Here, E12 = E„i/A is effective interbit-coupling with Em = 
4e^Cm/Cs. The effective charge energy E'^^ is iS^'' = 
Ec, (-1/2 + ng^)+Em (-1/4 + ngj2) (j, fc = 1, 2) with 
"9. = iCg.Vg^ + Cp^Vp)/{2e) and Ec, = ie^C^JC^. 
The effective Josephson energy of the SQUID is E j — 
2e^P cos{T:^j/^o) with Josephson energy Ej-'-' of the sin- 
gle junction. Above, (7s = CsjCsg — C^, and C^. is the 
sum of all capacitances connected to the jth box. The pesu- 
dospin operators are defined as ai"*^ — \0j) {Oj | — | Ij) (Ij | and 
a« = |0,)(l,| + |l,)(0,|. 

First, let us consider the circuit working at the co-resonance 
point (i.e., Ug^ = Ug^ = 0.5, yielding E^^ — E^^ — 0), 
and the applied fluxes are set as ^j=0, k ^ j 

(yielding Ej^^—2ej \ Ej^^'—O). In this case, the circuit has 
the Hamiltonian 

i/i = -4^V« + £;i2a«af . (2) 

The corresponding time-evolution operator reads 



ergy. The two Cooper-pair boxes are coupled via the ca- 
pacitance Cm. The qubits work in the charge regime, with 



fesT < E'^^ < E^(^' < A(j = 1,2), where both quasi- 
particle tunnelling and excitations are effectively suppressed 
and the number rij (with nj — 0, 1, 2, ...) of Cooper-pairs in 

the boxes is a good quantum number. Here, ks, T, A, E^\ 

and Ej are the Boltzmann constant, temperature, supercon- 
ducting gap, the charging and Josephson energies of the jth 
qubit, respectively. Following Refs. [4j(5J, the system is as- 
sumed to work near the co-resonant point and its quantum 
dynamics can be restricted to the subspace spanned by the 
four lowest charge states: |00), |10), |01) and |11). Thus, the 
Hamiltonian of this circuit is 



H 



-El 



E12 C7(i)a(2) 



(1) 



Ui {t) = exp 



it ^ 



exp 



Uintit), (3) 



where the operator C/int {t) is determined by 



ih 



dUintjt) 

dt 



= HnJt)Uint{t), 



with 



Hint{t) = -Biaexp 



Considering Q = Ei2/{'ie[p) < 1 (e.g., Cj < 1/4 in the cir- 
cuit t5J), one can make the following perturbation expansion: 



ai^Vf)exp 



(4) 



f/intW = 1 + 
= 1- 



dt'H„At') 



2 „t rt' 



dt'dt"Hi^t{t')Hint{t") 



2e 



(5) 



Neglecting the higher-order terms of Q, since it is small, the 
Hamiltonian of the system can be effectively rewritten as 



H 



off 



^12 



2e 



J J 



(6) 



the time evolution, due to its small probability (proportional 
to Q). Under this approximation, the fixed interaction 
between the qubits has been effectively eliminated, except 
resulting in a shift of the relatively strong Josephson energy. 
Thus, the system effectively undergoes an evolution: 



Above, the first-order expansion term ul^l{t) = 
{—i/h)j'^dt'Hint{t') practically does not contribute to 



exp 



K 



(7) 
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FIG. 2: Probabilities Pappr (solid lines) and Pcx (dotted lines) for the 
transition lljOfc) ^ |OjOfe) versus r = e^j^t/h, due to the approxi- 
mate and exact time evolutions, respectively. The difference between 
them decreases with decreasing interbit coupling Q. 



which reduces to the single-qubit cri -rotation (i.e. , qubit-flip) 
on the jth qubit, if the duration is set by cos (pj = 0. 

The robustness of this dynamical decoupling can be verified 
by testing the difference of the corresponding physical effects, 
e.g., the transition probabilities P between two selected states, 
due to the present approximate time-evolution 



(8) 



and the exact one 

U^^{t) = exp(-zii/i/?i) = A(t)ai^') ® /C^) 
+ S(t)(|0,Ofc)(0,Ofc| + |l,lfe)(l,lfe|) 
+ i3*(i)(|l,0fc)(l,0fc| + |0,U)(0,l,|), 



(9) 



respectively. Here, 

A{t) = tp,{t), B{t) - [1 - p2(t)]i/2 exp[-zC,(t)], (10) 
with 

Pji^) = vj^ sm{e[pvjt/h), 

ij(t) = arctan[2Cj^'7nan(£yVji/ri)]. 

Figure 2 shows that the probabilities for the transition 
|ljOfc) |OjOfc) due to the evolutions C/appr(i) (solid lines) 
and Ucx{t) (dotted lines) possess the same oscillating period. 
Also, the difference between the probabilities decreases when 
deceasing the coupling strength Q : the largest differences are 
less than 0.06 and 0.04 for coupling strengths Q — 1/8 and 
Q = 1/10, respectively. 

Similarly, if the system works far from the co-resonance 



to the charging energy E'^'' . Thus, the Hamiltonian (1) can be 
effectively approximated by 



(11) 



With £;,=4^')[l + c;|/(l-q22)], and ^,-i?yV(24'')' ^12 = 
Ei2/Ej^\ The evolution corresponding to this effective 
Hamiltonian results in a two-qubit operation 



i?(12)(;^)^CXp[-ZXl2a(lVp)] 



5 Jl exp[-tXj(T 

J = l,2 



(12) 



where x = {Xj^Xu} and Xj 
For the simplest case where E^j 

% ■ 



(i) 



Ejt/h, xi2 = Ei2t/h. 
— and thus c,j — 0, 



(7) 

we have Ej = E^ . By using a refocusing technique, like 



in NMR fy, we can effectively realize another important 
single-qubit operation: 



(13) 



with — 2xj- The inverse of this operation, i.e., the gate 

Ri^\^4'j) = exp [if/jjcr can be obtained by changing the 
signs of Ej via controlling the applied gate voltage. 

The single-qubit gates Rx ^ (ipj) and {(f>j) do not com- 
mute, and thus constitute a universal single-qubit gate set, 
which can assist the realization of two-qubit gates to imple- 
ment any unitary operation on this circuit. For example, a 
Hadamard-like operation 



R,{e,) = i?(^-)(0,/2)4^)(7r/4)4^)(-0,/2) 
_ 1 / 1 —iexp{i6j) 



1 



(14) 



can be implemented, which will take an important role for 
testing Bell's inequality. 



III. TESTING BELL'S INEQUALITY BY USING 
EFFECTIVE SINGLE-QUBIT LOCAL OPERATIONS 

By using the above dynamical decoupling procedure, we 
now show that Bell's inequality may be tested with fixed- 
coupling Josephson-qubits. 

First, the desired entanglement between these SQUID- 
based qubits can be created in a repeatable way. Initially, 
the system works sufficiently far from the co-resonance point 
and remains at the state \ipo) — |00), $j = 0. Now, a pair 
of gate voltage pulses brings the system to the co-resonance 
point 1 5] and lets the system undergo the evolution U3{t) — 
exp{—itH3/h), with 

^3 = ~ej J2 ^i'^ + E12 <J^P ® of^ . (15) 

J = l,2 



point (e.g., Ug. < 0.25, |4]), then both E j and the fixed For simplicity, here we assume that e j 



(1) 



-(2) 



ej. We 



coupling E12 are the small perturbative quantities, compared analytically derive the time-dependent degree of entanglement 
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or concurrence f?! Ce (t) of this circuit 



with 



P{t) = cos^ ^{t) - cos g{t) 



sin^ d{t) 



(16) 



and 



Q{t) 



m 
lit) 



sin g{t), 



(17) 
(18) 



1 + c- 



7(t)(i + C')'/', git) = 2C7W, 

2£jV^, C - E,2/i2ej). 



Figure 3 shows this evolution, showing some plateaus near the 
times te when sini?(te) — 0. At these times, the system is in 
the following compact entangled state 



with 



and 



m=am+(3\U), (19) 
CEite) = 2\a+a_\ = \sm{E,2te/fi)l (20) 



a± = [l±exp{±it,Ei2/ti)]/2 , 

These states are very adjacent to the eigenstates of H^, and 
almost do not evolve for several very short time intervals. 
Thus, periodically, their degrees of entanglement are almost 
unchanged, shown by the short plateaus in Fig. 3. The max- 
imally entangled states (corresponding to the top plateaus in 
Fig. 3) occur when the pulse durations are set properly such 
that the condition coa{Ei2te/fi) = is further satisfied. 

Next, using dynamically-generated single-qubit operations, 
the controllable variables {6j} can be encoded into the gen- 
erated entangled states, keeping the degree of entanglement 
unchanged. The change of concurrence of the two-qubit en- 
tangled state can be effectively suppressed by continuously 
applying controllable single-qubit operations. This is similar 
to the approaches for suppressing decoherence in open quan- 
tum systems by using the quantum Zeno effect and the "bang- 
bang" decoupling method L16.I . Thus, the new entangled state 



has the same degree of entanglement as that of the \ipe) gen- 
erated above. 

Finally, the correlations between the classical variables 
{6j} can be measured by simultaneously detecting the pop- 
ulations of qubits in the excited |1) or ground states |0) 1 17ll. 
Experimentally, the above steps can be repeated many times 
for the same classical variables and then the correlation func- 
tion Ee{9i,62) can be measured as 



E,, 



iVsame(gl,g2)-iVdiff(gl,g2) 
Ns^meiei,02)+N^isieu02)' 



(22) 



with A^samc (^^1 , 02 ) (Ndis i^i ,02)) being the number of events 
with two qubits found in the same (different) logic states. The- 
oretically, the above projected measurements can be expressed 
via 

Pt = |11>(11| + |00)(00| - |10)(10| - |01)(01| 

a'i^®af\ (23) 

and the correlation in the outcomes can be calculated as 



E{0,,02) - mPrWe) 

— ± sin(fei?i2//i) sin(6'i 



(24) 



For the sets of angles: {0^,0'^} = {~2l/^' 37r/8}, the 
Clauser, Horne, Shimony and Holt (CHSH) @] function 

/(IV-^)) = \E{0,,02) + E{0[,e2) + E{e,,e'2)-E{e[,0'^)\ 

= 2V2 \sm{teEi2/h)\ (25) 



is larger than 2 for 



\sin{t,Ei2/h)\ > 



V2- 



(26) 



Therefore, properly setting the pulse duration to prepare 
the desired entangled state (whose plateau-like concurrence is 
lai-ger than 0.707), the CHSH-type Bell's inequahty 0] 



feiWe)) < 2 



(27) 



IV'e) = n RM)\^^ 



can be effectively tested by experimentally measuring the 
CHSH function: fei^) = \Ee{0i,02) + E,{0[,62) + 

Eei01,0'2) — Ee{0'i,92)\- 



IV. DISCUSSIONS AND CONCLUSIONS 



The simplest dynamical-decoupling process proposed 



rn.n— 0,1 



1 1 mn) , (21) above consists of two ax -pulses and two delays 



with 



Udir) = exp[-iEi2Tal^>a 



(1)^(2)1 



(28) 



aoo = [a - f3exp{i0i + i62)]/2, 

aio = [—iaexp{~i9i) ~ i/S exp{i02)]/2, 

ttQi = [—iaexp{—i92)~i(3exp{i0i)]/2, 

ail = [aexp{-i0i - i02) - l3]/2, 



The duration of the tri'' ■'-pulse is calculated fs*] as tx={2l + 

l)to, l^0,l,2,...;andto=TTh/[2e^p{l + 2Cf)] ~ 31 ps. Con- 
sequently, the longest time delay t between the two cTj; -pulses 
could be estimated as r 270 ps (for 1^0), - 200 ps (for 
1=1), - 145 ps (for 1=2), and - 98 ps (for 1=3), etc. Thus, 
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FIG. 3: Dynamical evolution of the degree of entanglement Csit) 
for the circuit (Fig. 1) with couplings: Em/£j ~ 1/2 (dashed 
line), 1/4 (solid line), respectively. Here, ej — 30 /leV. The 
Ce = 0.707 (dotted line) gives the threshold for the violation of 
Bell's inequality for the entangled state (8), whose degree of entan- 
glement slowly changes during several short time intervals (see the 
plateaus in the figure). 

it is easy to experimentally check this simplest proposal for 
eliminating the fixed interbit coupling by using the pulse se- 
quence: UdiT)ai'''Ud{T)ay\ OT ay^UdiT)ay^Lfd{T). After 
these operations the two qubits should return to their initial 
states. Furthermore, an universal two-qubit controlled-CTz 
gate could be implemented by using the operational sequence: 
Ud{-7T/m2)Ri''\7T/A)Ri'\n/A). 

Similar to other theoretical schemes (see, e.g., 1 18|), the re- 
alizability of the present proposal also faces the technological 
challenge of rapidly on/off switching the Josephson energy of 
the qubit by using fast magnetic pulses 1 19]. This experimen- 
tal difficulty could be relaxed by increasing the durations of 



the applied pulses. Especially, the decoherence time of the 
two-qubit capacitively-coupled Josephson circuit |20] could 
be increased by decreasing the coupling capacitance Cm- In 
principle, the lifetime of the generated entangled state (19) ad- 
equately allows to perform the required operations for testing 
Bell's inequality (27), since such a state is very adjacent to 
the eigenstates of the circuit's Hamiltonian ^^3. In fact, the 
decay time of a two-qubit excited state is long (up to ~ 0.6 

(i) 

ns), even for very strong interbit coupling (e.g., Em — '2s j 
in the recent experiment fs*]). In addition, the influence of 
the environmental noise and operational imperfections is not 
fatal, as the nonlocal correlation E{9i,6j) in Bell's equality 
is statistical — i.e., its fluctuations could be effectively sup- 
pressed by averaging over several repeatable experiments. 

In summary, we propose an effective dynamical-decoupling 
approach to overcome the fixed-interaction difficulty in super- 
conducting nanocircuits. The dynamically-generated single- 
qubit operations may be used to test Bell's inequality, pro- 
viding another way to verify the existence of entanglement 
between two capacitively coupled Josephson qubits. The pro- 
posed approach can be easily modified to manipulate quantum 
entanglement in other "fixed-interaction" solid-state systems, 
e.g., the capacitively (inductively) coupled Josephson phase 
(flux) circuits, and the Ising (Heisenberg)-spin chain. 
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